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The 3-D Z(3) Potts model is a model for finite temperature QCD with heavy quarks. The chem- 
ical potential in QCD becomes an external magnetic field in the Potts model. Following Alford 
et al.[g|], we revisit this mapping, and determine the phase diagram for an arbitrary chemical po- 
tential, real or imaginary. Analytic continuation of the phase transition line between real and 
imaginary chemical potential can be tested with precision. Our results show that the chemical 
potential weakens the heavy-quark deconfinement transition in QCD. 
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1. Introduction 

The 3 -dimensional 3-state Potts model and the finite temperature QCD with infinitely heavy 
quarks share the same Z(3) global symmetry, and in both theories the order parameter shows a 
first order phase transition. When an external magnetic field is turned on in the Potts model, this 
Z(3) symmetry is explicitly broken and the first order phase transition is weakened, becoming 
second order for a certain strength of the magnetic field and crossover beyond. In the case of finite 
temperature QCD, the quark mass plays the role of the external magnetic field in the 3-D 3-state 
Potts model. The first order phase transition of the quenched theory (infinitely heavy quark limit) 
weakens, and turns into a second order transition for a critical, heavy but not infinite, dynamical 
quark mass [|l|, ^J. A universality argument then implies that the critical properties of the two 
theories, Potts and QCD, will be the same at this second-order point. Since the 3-D 3-state Potts 
model is simpler than QCD and can be simulated on a large lattice, its numerical investigation will 
give us more detailed quantitative information on the critical properties of finite temperature QCD 
with heavy quarks. 

The effect of a quark chemical potential on hot QCD with static quark sources can be formu- 
lated in the following grand canonical ensemble [^] : 

Z(ji) = ^Z nTl e^"^ 

n,n 

= V f [9U] — (<£>[U]) n ^- ($[[/]*)" e Sg[U]-Pn(M-n)-Pn(M+n) 

= J[@U}exp(-S f ,[U}+e- p{M -L l) <i>[U}+e- p{M+ ^<i>[U}*), (1.1) 

where Z„ ^ is the canonical partition function with n quarks and n anti-quarks, U is the SU(3) gauge 
field, S g is the gauge action, <J> is the Polyakov line, <!>* the anti-Polyakov line, M is the heavy 
quark mass and /I is the quark chemical potential. When there is no chemical potential (jj. = 0), the 



action in (|L1|) is real. If \x / 0, the action becomes complex and Monte Carlo simulation is difficult 
because of the "sign problem": the usual probabilistic interpretation of Eq. ( |1.1[ ) is not possible. 

Symmetry consideration tells us that the important dynamics of the gauge field are those of 
the Polyakov line and the anti-Polyakov line. Thus, the corresponding lattice Hamiltonian for the 
3-D 3-state Potts model is given by: 

H = -k^8^ )M , +i) - Y,[ha>(fi+hfQ*W] (1.2) 

i,x x 

where h = e -^ M -^ = h M e^ = h M e* and ti = e -/»(*+M) = h u e~^ = h M e^ with h M = e- M l T , 
/T = n/T, and $(x) is a Z(3) spin. If h ^ h'*, this Hamiltonian also is complex. However, the 
partition function remains real [ph. The partition function includes summation over all the possible 
Z(3) spin configurations. Introducing "bonds" with a certain probability among parallel Z(3) spins 
and defining a "cluster" made of the sites connected by the bonds, we can divide the summation 
into the sum over clusters and that within a cluster. Then, after analytically summing over Z(3) spin 
orientations within each cluster, the partition function in terms of cluster configurations is given by 



Z= j \®b\(e k - l f"Y[ p / ' M|C|coshM + 2e- ftM|c|cosh ' I cos(y3/jM|C|sinhAi) 



(1.3) 
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where 3>b is the sum over all the possible bond configurations, Nb is the number of bonds in 
a given cluster configuration, and \C\ is the number of sites in a given cluster. Z is real and is 
free from the "sign problem" since the second term which can be negative due to the presence of 
cos(\/3h M \C\ sinh/T) is always smaller than the first term. However, the "solution" of the complex 
action problem in the 3-D 3-state Potts model is different from the SU (2) gauge theory case [||] 
and from the four-fermi theory case [0], where the action itself can be shown to be real even with a 
real chemical potential. In the Potts model, the action is complex, and a change of variables (from 
spins to bonds) is necessary to recover a real effective action and show that the partition function 
of the model remains real. 

The (h = h!) case (that is, zero chemical potential) has been studied in [Q]. The (h' = 0) case 
(that is, M,n — » oo while Mj\i is fixed) has been investigated in [Hp. Since we are also interested in 
testing the analytic continuation of imaginary chemical potential results to real chemical potential, 
instead of taking h' = limit, we study here the full parameter space of arbitrary (h,h'). 

For the case of an imaginary chemical potential {}X — ► i.e., h! = liMe~ liil = h*), the partition 
function, Zj, becomes 



f 9b{e k - l) Nb Y[ \e 2flM ^ C0S ^ +2e- /lM l c l cos ^cosh(V3/j M |C|sin7i / ) 



(1.4) 



Since cosh(i/i) = cos(/i),sinh(//i) = /sin(/i) and cos(/0) = cosh(0), the relation Z(jj, —> ijXi) = 
Zi(hi) is quite obvious and the analytic continuation between Z and Zj is transparent. The Roberge- 



Weiss symmetry Z/(^r) = Zj{^r + ir w )> can a l so be seen clearly in (1.4). 
2. Real Chemical Potential 



Using the Swendsen-Wang cluster algorithm []8J] on (|L3j), we simulate the 3-D 3-state Potts 
model. The actual simulation is performed along the h = h! line in the (h,h') parameter space and 
the data are re- weighted for arbitrary (h,h'). Although the reweighting factor is not always positive, 



this strategy is more efficient than sampling Eq. ( |1.3[ ), because the ensuing sign problem is very 
mild. Simulations are performed on lattice volumes 56 3 ,64 3 , and 72 3 . Typically ~ 2 million data 
sample is collected. The critical point is located by requiring the third order Binder cumulant of 
the Z(3) spin magnetization to vanish and by evaluating the fourth order Binder cumulant (B4) 
of the magnetization at that point. Since the universality class of the 3-D 3-state Potts model 
is that of the 3-D Ising model, B4 for the Potts model at the critical point should be equal to 
that of the Ising model (= 1.604(1)). Figure ([[]) shows the fourth order Binder cumulant as a 
function of h (here, h' = h). From this Figure, we find that the critical end point is at (k c ,h c ) = 
(0.54940(4)0.000255(5)) in comparison with the value given in [@|, (0.54938(2), 0.000258(3)). 
Similarly, for the h' = case, we obtain that (k c ,h c ) = (0.54947(1), 0.000465 (5)) ([§] reported 
(k c ,h c ) = (0.549463(13), 0.000470(2)). 

Similarly, the critical end point can be located for arbitrary h and h! . Figure (||) shows the Qi,h') 
parameters for the second order phase transition. Here, instead of directly plotting (h,h'), we use 
variables M/T and jl/T, which can be related to QCD. The line in the figure is from the asymptotic 
value, (h = h c ,h' = 0). Since h = e~ M l T+ti / T , the relation h c = e~ M / r +^/ r for given h c defines a 
line in the parameter space (y, f). In Figure (3), the simulation data lies above this asymptotic 
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Figure 1: B4 for magnetization vs h (h = h!) Figure 2: B4 for magnetization vs h (h 1 = 0) 



line. This means that the anti-Polyakov line in finite temperature QCD plays an important role in 
critical parameter space and should not be neglected in considering the critical properties of the 
theory. To further compare our results with the limit studied in we show T/M vs. jl/M in 
Figure (4). As jti/M increases, our simulation data approach the asymptotic line as expected. This 
figure clearly shows that the range of T/M values for which the transition is first-order shrinks as 
the chemical potential is turned on. In fact, the transition disappears altogether for ^ > 1. So the 
effect of the chemical potential is to weaken the phase transition. 



Potts, 72 A 3 Potts, 72 A 3 




(fl/T) A 2 H/M 

Figure 3: M/T for 2nd order transition vs. (jx/T) 2 Figure 4: T/M for 2nd order transition vs. jx/M 



Figure <J3]> suggests an interesting comparison with the finite temperature QCD phase diagram. 
Let us imagine increasing the chemical potential from zero, following the arrow in Figure (5). 
Without chemical potential, the quark mass is chosen such that the transition is first order. With 
a non-zero chemical potential, this first order phase transition weakens. At a certain chemical 
potential, the system shows a second order phase transition. With chemical potential larger than 
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this critical value, the system shows a cross-over. On the other hand, consider Figure (6) which is 
analogous to Figure (5). This figure is a schematic phase diagram for three-flavor QCD with light 
quarks suggested in ^ and summarizes the conventional wisdom. Without chemical potential, 
3-flavor QCD shows a first order phase transition when the quark mass is smaller than a certain 
mass (m c (0)). Let us say that we start with a system where the quark mass is larger than m c (0). 
Then, this 3-flavor QCD system has a cross-over. In this case, increasing the chemical potential 
makes the transition stronger: a second order phase transition appears when the chemical potential 
hits a certain magnitude. If the chemical potential is increased further, the system undergoes a first 
order phase transition. Therefore, the 3 light flavor QCD system shows (cross-over — ► second order 
phase transition — > first order phase transition) as the chemical potential is increased. In contrast, 
the heavy quark QCD system suggested by 3-D Z(3) Potts model will undergo (first order phase 
transition — ► second order phase transition — > cross-over) as the chemical potential is increased. 
Interestingly, some recent light-quark QCD simulations give support to the possibility that the 



conventional wisdom scenario above might be at fault Q10|], so that the influence of the chemical 
potential in the light-quark and the heavy-quark cases would be similar. 
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Figure 5: M/T for critical point vs. (j) 2 



Figure 6: schematic QCD phase diagram for 3 light quarks 



3. Imaginary vs Real Chemical Potential 



For the case of an imaginary chemical potential, one can perform a direct sampling of Eq. (|1.2|), 
or reweight the h = h' data used for the real chemical potential study. As before, the critical point 
is located by use of the magnetization Binder cumulant and B4 = 1.604 crossing point. 

In Figure (||), we put together the critical end point parameters obtained from the imaginary 
chemical potential case and those from the real chemical potential case (note the similarity with the 
schematic phase diagram for 3 light flavor QCD shown in Figure (6)). Although the proximity of 
the Roberge-Weiss transition at jit/ = n/3 introduces curvature in the imaginary chemical potential 
result, the small real /I result is smoothly connected to the small imaginary result and the analytic 
continuation poses no serious problem even on a large lattice volume such as 72 3 . However, the crit- 
ical line shows significant curvature, which limits the accuracy of the extrapolation from imaginary 



166/5 



The 3-states Potts model 



Seyong Kim 



to real p. A 4th order polynomial f = 8.273 + 0.585(£) 2 - 0.174(£) 4 + 0.160(£) 6 - 0.071(£) 8 
is necessary to describe the critical line from (j) 2 = —(f) 2 to 1.5. 

4. Summary 

We extend earlier work on the 3-D Z(3) Potts model with one external field coupled to <S>(3c) + 
<t>*(x) [Q] or to <S>(3c) []3|] to the full general case M>(x) +h'<t>*(x), and present a complete picture 
of the phase diagram in the (K,h,h') parameter space. Our investigation also gives us a handle 
on heavy quark QCD at finite density and temperature which shares the same global symmetry, 
allowing us to make statements about the heavy quark QCD phase diagram in the finite temperature, 
finite density, and heavy quark mass parameter space. 

For a real chemical potential, the partition function of the Potts model is shown to be real even 
though the action itself is complex. The sign problem is mild, and simulation results show that 
turning on a chemical potential pL makes the finite temperature transition weaker, so that the region 
of parameter space corresponding to a first-order transition shrinks under the influence of /J.. This 
implies that a similar phenomenon occurs in the phase diagram of heavy quark QCD. 

For an imaginary chemical potential, the action of the Potts model is real and the model is easy 
to study. Since both real chemical potential and imaginary chemical potential can be simulated, 
analytic continuation from imaginary to real chemical potential can be tested. We find that analytic 
continuation works satisfactorily, even with large lattice volumes such as 72 . 

In short, the 3-D Z(3) Potts model has a rich structure and provides us with a useful "proving 
ground" for studying the finite temperature/density phase structure of QCD. 
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